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this manifold in terms of its geodesic radius. We give several matching expressions 
for this fundamental solution including a definite integral over reciprocal powers of 
the hyperbolic sine, finite summation expression over hyperbolic functions, Gauss 
hypcrgcomctric functions, and in terms of the associated Lcgcndre function of the 
second kind with order and degree given by d/2 — 1 with real argument greater than 
unity. We also demonstrate uniqueness for a fundamental solution of Laplace's equation 
on this manifold in terms of a vanishing decay at infinity. 
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1. Introduction 

We compute closed-form expressions of a spherically symmetric Green's function 
(fundamental solution) of the Laplacian (Laplace-Beltrami operator) on a Riemannian 
manifold of constant negative sectional curvature, namely the hyperboloid model of 
hyperbolic geometry. Useful background material relevant for this paper can be found 
in Vilenkin (1968) [25], Thurston (1997) [23], Lee (1997) [18] and Pogosyan & Winternitz 
(2002) [22]. 

This paper is organized as follows. In section [2] we describe the hyperboloid model 
of hyperbolic geometry and its corresponding metric, global geodesic distance function, 
Laplacian and geodesic polar coordinate systems which parametrize points in this model. 
In section [3] for the hyperboloid model of hyperbolic geometry, we show how to compute 
radial harmonics in a geodesic polar coordinate system and derive several alternative 
expressions for a radial fundamental solution of the Laplacian on the d- dimensional R- 
radius hyperboloid with R > 0. In section H] we prove that our derived fundamental 
solution is unique in terms of a vanishing decay at infinity. 

Throughout this paper we rely on the following definitions. For 01,02, ... G C, 
if i,j G Z and j < i then J2n=i a n = an d 111=1 a " = 1- The set of natural 
numbers is given by N := {1, 2, . . .}, the set N := {0, 1, 2, . . .} = N U {0}, and the set 
Z := {0, ±1, ±2, . . .}. The set R represents the real numbers. 

2. The hyperboloid model of hyperbolic geometry 

Hyperbolic space in (^-dimensions is a fundamental example of a space exhibiting 
hyperbolic geometry. It was developed independently by Lobachevsky and Bolyai 
around 1830 (see Trudeau (1987) [21]). It is a geometry analogous to Euclidean 
geometry, but such that Euclid's parallel postulate is no longer assumed to hold. 

There are several models of <i-dimensional hyperbolic geometry including the Klein 
(see Figured]), Poincare (see Figure |2]), hyperboloid, upper-half space and hemisphere 
models (see Thurston (1997) [23]). The hyperboloid model for ^-dimensional hyperbolic 
space is closely related to the Klein and Poincare models: each can be obtained 
projectively from the others. The upper-half space and hemisphere models can be 
obtained from one another by inversions with the Poincare model (see section 2.2 in 
Thurston (1997) [23]). The model we will be focusing on in this paper is the hyperboloid 
model. 

The hyperboloid model, also known as the Minkowski or Lorentz models, are models 
of d-dimensional hyperbolic geometry in which points are represented by the upper sheet 
(submanifold) S + of a two-sheeted hyperboloid embedded in the Minkowski space R^ 1 . 
Minkowski space is a (d + l)-dimensional pseudo- Riemannian manifold which is a real 
finite-dimensional vector space, with coordinates given by x = (xq, x\, ■ ■ ■ , Xd)- It is 
equipped with a nondegenerate, symmetric bilinear form, the Minkowski bilinear form 

[x, y] = x y - x x yi - ... - x d y d . 
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Figure 1. This figure is a graphical depiction of stereographic projection from the 
hypcrboloid model to the Klein model of hyperbolic space. 

The above bilinear form is symmetric, but not positive-definite, so it is not an inner 
product. It is defined analogously with the Euclidean inner product for R d+1 

(x, y) = x y + xiyi + • • • + x d y d . 

The variety [x, x] = x\ — x\ — . . . — x 2 d = R 2 , for x £ R^' 1 , using the language of 
Beltrami (1869) [3] (see also p. 504 in Vilenkin (1968) [25]), defines a pseudo-sphere of 
radius R. Points on the pseudo-sphere with zero radius coincide with a cone. Points on 
the pseudo-sphere with radius greater than zero lie within this cone, and points on the 
pseudo-sphere with purely imaginary radius lie outside the cone. 

For R £ (0, oo), we refer to the variety [x, x] = R 2 as the i?-radius hyperboloid H^. 
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Figure 2. This figure is a graphical depiction of stercographic projection from the 
hypcrboloid model to the Poincarc model of hyperbolic space. 



This variety is a maximally symmetric, simply connected, d- dimensional Riemannian 
manifold with negative-constant sectional curvature (given by —1/R 2 , see for instance 
p. 148 in Lee (1997) [IE]), whereas Euclidean space R d equipped with the Pythagorean 
norm, is a space with zero sectional curvature. For a fixed R G (0,oo), the i?-radius 
hypersphere is an example of a space (submanifold) with positive constant sectional 
curvature (given by 1/R 2 ). We denote unit radius hyperboloid by H d := Hf and the 
unit radius hypersphere by S d := Sf. 

In our discussion of a fundamental solution for the Laplacian in the hyperboloid 
model of hyperbolic geometry, we focus on the positive radius pseudo-sphere which can 
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be parametrized through subgroup-type coordinates, i.e. those which correspond to a 
maximal subgroup chain 0(d, 1) D ... (see for instance Pogosyan & Winternitz (2002) 
|22j). There exist separable coordinate systems which parametrize points on the positive 
radius pseudo-sphere (i.e. such as those which are analogous to parabolic coordinates, 
etc.) which can not be constructed using maximal subgroup chains (we will no longer 
discuss these). 

Geodesic polar coordinates are coordinates which correspond to the maximal 
subgroup chain given by 0(d, 1) D 0(d) D .... What we will refer to as standard 
geodesic polar coordinates correspond to the subgroup chain given by 0(d,l) D 0(d) D 
0(d — 1) D ■•• D 0(2). Standard geodesic polar coordinates (see Olevskii (1950) 
[20J; Grosche, Pogosyan & Sissakian (1997) [13]), similar to standard hyperspherical 
coordinates in Euclidean space, can be given by 

xo = R cosh r 
Xi = R sinh r cos 9\ 
X2 = R sinh r sin 9\ cos 62 



Xd-2 = R sinh r sin 9 1 ■ ■ ■ cos Qd-2 
Xd-i = R sinh r sin Q\ ■ ■ ■ sin Qd-2 cos 1 
Xd = R sinh r sin #1 • • • sin Qd-2 sin < 



(1) 



d-2}. 



where r G [0, 00), £ [0, 27r), and 8 { G [0, 7r] for % G {1, 

The isometry group of the space is the pseudo-orthogonal group SO(d, 1), the 
Lorentz group in (d + l)-dimensions. Hyperbolic space H^, can be identified with the 
quotient space SO(d, l)/SO(d). The isometry group acts transitively on H^. That is, 
any point on the hyperboloid can be carried, with the help of a Euclidean rotation of 
SO(d — 1), to the point (cosh a, sinh a, 0, . . . , 0), and a hyperbolic rotation 

= —x\ sinh a + xq cosh a 
- —x\ cosh a — Xq sinh a 



Xr 



L 1 



maps that point to the origin (1,0,..., 0) of the space. In order to study a fundamental 
solution of Laplace's equation on the hyperboloid, we need to describe how one computes 
distances in this space. 

One may naturally compare distances on the positive radius pseudo-sphere through 
analogy with the i?-radius hypersphere. Distances on the hypersphere are simply given 
by arc lengths, angles between two arbitrary vectors, from the origin, in the ambient 
Euclidean space. We consider the d- dimensional hypersphere embedded in R d+1 . Points 
on the hypersphere can be parametrized using hyperspherical coordinate systems. Any 
parametrization of the hypersphere S^, must have (x, x) = Xq + . • . + x 2 d = R 2 , with 



R > 0. The distance between two points on the hypersphere x, x' G is given by 



d(x, x') = i?7 = R cos" 



X, X 



:x, x )(x',x'; 



Rcos 1 ( — (x,x'; 



(2) 
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This is evident from the fact that the geodesies on are great circles (i.e. intersections 
of with planes through the origin) with constant speed parametrizations (see p. 82 
in Lee (1997) [18]). 

Accordingly, we now look at the geodesic distance function on the <i-dimensional 
positive radius pseudo-sphere H^. Distances between two points on the positive radius 
pseudo-sphere are given by the hyperangle between two arbitrary vectors, from the 
origin, in the ambient Minkowski space. Any parametrization of the hyperboloid H^, 
must have [x, x] = R 2 . The geodesic distance between two points x, x' £ is given 
by 

4x,x<) = ficosh- ygjj) = Hcosh- (i,[x,xl) , (3) 

where the inverse hyperbolic cosine with argument x £ (1, oo) is given by (see (4.37.19) 
in Olver et al. (2010) [21]) 

cosh -1 x = log (^x + V x 2 — lj . 

Geodesies on are great hyperbolas (i.e. intersections of H^. with planes through the 
origin) with constant speed parametrizations (see p. 84 in Lee (1997) [H]). We also 
define a global function p : H d x H d — >• [0, oo) which represents the projection of global 
geodesic distance function <^fy on onto the corresponding unit radius hyperboloid 
H d , namely 

p(x,x') :=d(x,x')AR, (4) 
where x = x/i? and x' = x'/i?. 



2.1. The Laplacian on the hyperboloid model 

Parametrizations of a submanifold embedded in either a Euclidean or Minkowski space 
is given in terms of coordinate systems whose coordinates are curvilinear. These 
are coordinates based on some transformation that converts the standard Cartesian 
coordinates in the ambient space to a coordinate system with the same number of 
coordinates as the dimension of the submanifold in which the coordinate lines are curved. 

On a (i-dimensional Riemannian manifold M (a manifold together with a 
Riemannian metric g), the Laplace-Beltrami operator (Laplacian) A : C P (M) — > 
C P ~ 2 (M), p > 2, in curvilinear coordinates £ = . . . , C, d ) is given by 




where \g\ = \ det(gij)\, the infinitesimal distance is given by 

d 

ds 2 = J2 (hidCd?.. (6) 



Fundamental solution of Laplace 's equation in hyperbolic geometry 



7 



and 

d 



i=i 



where b~\ G {0, 1} is the Kronecker delta defined for all j,j G Z such that 
f 1 if % = ?', 

[0 l£ l y£ J. 

For a Riemannian submanifold, the relation between the metric tensor in the ambient 
space and of (0) and (J5]) is 

dx k dx l 



Hi AO = Yl Gk 



k,l=0 s s 

On the ambient space is Minkowski, and therefore Gij = diag(l, — 1, . . . , — 1). 

The set of all geodesic polar coordinate systems on the hyperboloid correspond 
to the many ways one can put coordinates on a hyperbolic hypersphere, i.e., the 
Riemannian submanifold U C defined for a fixed x' G such that rf(x, x') = b — 
const, where b G (0, oo). These are coordinate systems which correspond to subgroup 
chains starting with 0(d, 1) D 0(d) D • • •, with standard geodesic polar coordinates 
given by (PQ) being only one of them. (For a thorough description of these see section 
X.5 in Vilenkin (1968) [25].) They all share the property that they are described by 
(d + l)-variables: r G [0, oo) plus ci-angles each being given by the values [0, 2tt), [0, n], 
[-7r/2,7r/2] or [0,7r/2] (see Izmest'ev et al. (1999, 2001) HQ El). 

In any of the geodesic polar coordinate systems, the global geodesic distance 
between any two points on the hyperboloid is given by (cf. (j3J)) 

d(x, x') = R cosh -1 (cosh r cosh r' — sinh r sinh r' cos 7) , (8) 

where 7 is the unique separation angle given in each hyperspherical coordinate system. 
For instance, the separation angle in standard geodesic polar coordinates (QQ) is given 
by the formula 

d-2 d-2 i-1 

cos 7 = cos(0 — <p') Y[ sin ^sin 9/ + cos 0jCOS 9/ J^J sin 9jSm 9/. (9) 

i=l i=l j=l 

Corresponding separation angle formulae for any geodesic polar coordinate system can 
be computed using (J2]), ([3]), and the associated formulae for the appropriate inner- 
products. Note that by making use of the isometry group SO(d, 1) to map x' to the 
origin, then p = Rr for and in particular p = r for H d . Hence, for the unit 
radius hyperboloid, there is no distinction between the global geodesic distance and 
the r-parameter in a geodesic polar coordinate system. For the i?-radius hyperboloid, 
the only distinction between the global geodesic distance and the r-parameter is the 
multiplicative constant R. 
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The infinitesimal distance in a geodesic polar coordinate system on this submanifold 
is given by 

ds 2 = R 2 (dr 2 + sinh 2 r d 7 2 ), (10) 

where an appropriate expression for 7 in a curvilinear coordinate system is given. If one 
combines (PQ), (jSJ), © and (fTOl . then in a particular geodesic polar coordinate system, 
Laplace's equation on is given by 



d 2 f df 1 

+ (d - 1) cothr^- + , 2 A S d-i/ 



0, fir 



<9r 2 <9r sinh 2 r 

where Agd-i is the corresponding Laplace-Beltrami operator on the unit radius 
hypersphere S d_1 . 

3. A Green's function in the hyperboloid model 

3.1. Harmonics in geodesic polar coordinates 

Geodesic polar coordinate systems partition the i?-radius hyperboloid into a family 
of (d— l)-dimensional hyperbolic hyperspheres, each with a radius r G (0, 00), on which 
all possible hyper spherical coordinate systems for S^" 1 may be used (see for instance 
Vilenkin (1968) [25]). One then must also consider the limiting case for r = to fill out 
all of H^. In geodesic polar coordinates one can compute the normalized hyperspherical 
harmonics in this space by solving the Laplace equation using separation of variables 
which results in a general procedure which is given explicitly in Izmest'ev et al. (1999, 
2001) [TBI ITT] . These angular harmonics are given as general expressions involving 
trigonometric functions, Gegenbauer polynomials and Jacobi polynomials. 

The harmonics in geodesic polar coordinate systems are given in terms of a radial 
solution multiplied by the angular harmonics. The angular harmonics are eigenf unctions 
of the Laplace-Beltrami operator on S d ~ l with unit radius which satisfy the following 
eigenvalue problem 

A sd -tY l K (Z) = -l(l + d-2)Y l K (5l), 

where x G S^" 1 , Y t K (Si) are normalized hyperspherical harmonics, / G N is the angular 
momentum quantum number, and K stands for the set of (d — 2)-quantum numbers 
identifying degenerate harmonics for each I. The degeneracy 

(see (9.2.11) in Vilenkin (1968) [25]), tells you how many linearly independent solutions 
exist for a particular I value and dimension d. The hyperspherical harmonics are 
normalized such that 



S d-i 

where duo is the Riemannian (volume) measure (see for instance section 3.4 in Grigor'yan 
(2009) [12]) on S d_1 which is invariant under the isometry group SO(d) (cf. f )12p ). and for 
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x + iy = z G C, ~z = x — iy, represents complex conjugation. The generalized Kronecker 
delta 5ft (cf. (J7|)) is defined such that it equals 1 if all of the (d — 2)-quantum numbers 
identifying degenerate harmonics for each I coincide, and equals zero otherwise. 

Since the angular solutions (hyperspherical harmonics) are well-known (see Chapter 
IX in Vilenkin (1968) [25]; Chapter 11 in Erdelyi et al. (1981) [7]), we will now focus 
on the radial solutions on in geodesic polar coordinates, which satisfy the following 
ordinary differential equation (cf. (ITT]) ) for all R G (0, oo), namely 

d 2 u . , . , du 1(1 + d — 2) 

— + (d - 1) cothr- = -u = 0. 

dr 1 dr sinh r 

Four solutions to this ordinary differential equation u d ±, i*2± : (1; °°) ~~ C are given by 



and 



uJi(coshr) = . . d/2 _! Pdtt-l 1+0 (cosh r) 
sinh 1 r ' 



M2±(coshr) = . ud/2 _! Qd/2-l 1+0 (cosh r), 
sinh 1 r 



where Pjf,Q^ : (l,oo) — > C are associated Legendre functions of the first and second 
kind respectively (see for instance Chapter 14 in Olver et al. (2010) [2~T]). 

Due to the fact that the space is homogeneous with respect to its isometry 
group, the pseudo-orthogonal group SO(d,l), and therefore an isotropic manifold, we 
expect that there exist a fundamental solution of Laplace's equation on this space with 
spherically symmetric dependence. We specifically expect these solutions to be given 
in terms of associated Legendre functions of the second kind with argument given by 
coshr. This associated Legendre function naturally fits our requirements because it is 
singular at )r = and vanishes at infinity, whereas the associated Legendre functions 
of the first kind, with the same argument, are regular at r = and singular at infinity. 



3.2. Fundamental solution of the Laplacian 

In computing a fundamental solution of the Laplacian on H^, we know that 
-A?4(x,x') = 5 9 (x,x'), 

where g is the Riemannian metric on and <5 ff (x, x') is the Dirac delta function on 
the manifold H^. The Dirac delta function is defined for an open set U C with 
x, x' G such that 

f f 1 if x' G U, 

y^(x,x)dvoi 9 = j o . fx ,^ 

where dvo\ g is the Riemannian (volume) measure, invariant under the isometry group 
SO(d, 1) of the Riemannian manifold H^, given (in standard geodesic polar coordinates) 
by 

dvo\ g = R d sim/" 1 r dr du := R d sinh^ 1 r dr sm d ~ 2 B d _ x ■ • • sin 9 2 d6 l ■ ■ ■ d6 d ^. (12) 
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Notice that as r — > + that dvol s goes to the Euclidean measure, invariant under 
the Euclidean motion group E(d), in spherical coordinates. Therefore in spherical 
coordinates, we have the following 

S(r-r') S(9 1 -9' 1 )---5(9 d ^ 1 



5 3 (X,X 



'd-1, 



R d sinl/" 1 r> 



sin 6' 



'2 ■ ■ ■ sin 

In general since we can add any harmonic function to a fundamental solution of the 
Laplacian and still have a fundamental solution, we will use this freedom to make our 
fundamental solution as simple as possible. It is reasonable to expect that there exists a 
particular spherically symmetric fundamental solution H^(x, x') on the hyperboloid 
with pure radial p(x, x') := d(x, x')/i? (cf. (j4|)) and constant angular dependence 
(invariant under rotations centered about the origin), due to the influence of the point- 
like nature of the Dirac delta function. For a spherically symmetric solution to the 
Laplace equation, the corresponding Agd-i term vanishes since only the I = term 
survives. In other words, we expect there to exist a fundamental solution of Laplace's 
equation such that 7^(x, x') = f(p). 

We have proven that on the i?-radius hyperboloid H R , a Green's function for the 
Laplace operator (fundamental solution for the Laplacian) can be given as follows. 

Theorem 3.1. Let d E {2, 3, . . .}. Define X d : (0, oo) -> R as 

f°° dr 

Up) ■- 



d-2 qi 



sinh d 1 x ' 



(H* xH£)\{(x,x):xeH'} 

LWg) r ( ) 



R defined such that 



x, x' E H d R , and H R 

xo .- 27rd/2Rd 

where p := cosh -1 ([x, x']) is the geodesic distance between* and*! on the pseudo-sphere 
of unit radius H d , withx = x/R, x' — ~x!jR, thenl-i R is a fundamental solution for — A 
where A is the Laplace- Beltrami operator on H R . Moreover, 



Up) 



(d-2)!! 
(d 



d/2-l 

log coth - + cosh p 

2 (2* 



(2fc-2)! 



l)!!sinh 2fc 



P. 



if d even, 



\(d-l)/2 



+ 



3)!! 



(d-2)!! 

q \ (d~l)/2 



-lfcoth 



2fc— 1 



fc=l 



(2fc-l)(fc-l)!((d-2fc-l)/2)! 



> if d odd. 



or 



-l)/2 



(d-3)! 
(d-2)! 



(d-l)/2 

i + cosh p 

fc=l 



(2ife-3)!!( 



(2A;-2)!!sinh 2fe-1 p 
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^(ii, d+1 1 



(d — 1) cosh p smh d 2 p~ \2' 2 'cosh 2 p/' 
e -Mr(d/2-i) d/2 _ i 
= 2"/2-ir (d/2) sinh^ 1 p ^(coshp), 

where !! is i/ie double factorial, 2F1 is the Gauss hypergeometric function, and Q^, is the 
associated Legendre function of the second kind. 

In the rest of this section, we develop the material in order to prove this theorem. 

Due to the fact that the space H R is homogeneous with respect to its isometry 
group SO(d, 1), and therefore an isotropic manifold, without loss of generality, we 
are free to map the point x' G H R to the origin. In this case the global distance 
function p : H d x H d — > [0, 00) coincides with the radial parameter in geodesic polar 
coordinates, and we may interchange r with p accordingly (cf. (]§]) with r' = 0) in 
our representation of a fundamental solution of Laplace's equation on this manifold. 
Since a spherically symmetric choice for a fundamental solution of Laplace's equation is 
harmonic everywhere except at the origin, we may first set g = f in (fTTj) and solve the 
first-order equation 

g' + (d — 1) coth p g = 0, 

which is integrable and clearly has the general solution 

^(p) = ^ = c sinh 1 - d p, (13) 

where Cq G R is a constant which depends on d. Now we integrate (1131) to obtain a 
fundamental solution for the Laplacian in 

?4(x,x') = c I,(p) + Cl , (14) 



where 



p 



and Co, c\ G R are constants which depend on d. This definite integral result is mentioned 
in section II. 5 of Helgason (1984) [14] and as well in Losev (1986) [19] . Notice that 
we can add any harmonic function to ( TT4"1) and still have a fundamental solution of the 
Laplacian since a fundamental solution of the Laplacian must satisfy 



/ (-A^(x')?4(x,x0dvoi; = ^(x) 



for all cp G T>(H.j i ), where T> is the space of test functions, and cZvol^ is the Riemannian 
(volume) measure on H^, in the primed coordinates. In particular, we notice that from 
our definition of X d (TT5|) that 



lim X d (p) = 0. 
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Therefore it is convenient to set c± = leaving us with 

?4(x,x') = c X d (p). (16) 

In Euclidean space R d , a Green's function for Laplace's equation (fundamental 
solution for the Laplacian) is well-known and is given in the following theorem (see 
Folland (1976) |SJ ; p. 94, Gilbarg & Trudinger (1983) [9]; p. 17, Bers et at (1964) @], 
p. 211). 

Theorem 3.2. LetdeN. Define 

W2) , ll2 _ 



x'll^ if d= 1 or d> 3, 



— log ||x — x'H 1 if d = 2, 

27T 

i/ien Gr is a fundamental solution for —A in Euclidean space H d , where A is the Laplace 
operator in R d . 

Note most authors only present the above theorem for the case d > 2 but it is easily- 
verified to also be valid for the case d — 1 as well. 

The hyperboloid H^, being a manifold, must behave locally like Euclidean space 
R d . Therefore for small p we have e p ~ 1 + p and e~ p ~ 1 — p and in that limiting 
regime 

- log p if d = 2, 



pd-2 



if d > 3, 



which has exactly the same singularity as a Euclidean fundamental solution for Laplace's 
equation. Therefore the proportionality constant cq is obtained by matching locally to 
a Euclidean fundamental solution of Laplace's equation 

n% = c i d ~ g d , 

near the singularity located at x = x'. 

We have shown how to compute a fundamental solution of the Laplace-Beltrami 
operator on the hyperboloid in terms of an improper integral f|T5|) . We would now like to 
express this integral in terms of well-known special functions. A fundamental solution 
Id can be computed using elementary methods through its definition (|T5|) . In d = 2 we 
have 

J ( ) (°° ^ X — ^ 1 g COS ^~ ^ * — 1 th ^ 
7 „ sinh x 2 cosh p — 1 2 ' 



and in d = 3 we have 



sinh 2 x sinh p 



X s(p) = / — 2 = = cotn '° - L 
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This exactly matches up to that given by (3.27) in Hostler (1955) [T5]. In d G {4, 5, 6, 7} 
we have 

_ , . 1 , , p coshp 

X 4 (p)=--logcoth- + ^ ? -, 

Hp) = ^(coth 3 P -i)-(cothp-i), 

3 p cosh p 3 cosh p 

MP) = o lo § coth 9 + 7— 7a ~ . 7 2 , and 

8 Z 4 smh p 8 sinn p 

1 2 

J 7 (p) = -(coth 5 p - 1) (coth 3 p - 1) + cothp - 1. 

5 3 

Now we prove several equivalent finite summation expressions for X^(p). We wish 
to compute the antiderivative J m : (0, 00) —> R, which is defined as 

3m{x) := /shih^v 

where m G N. This antiderivative satisfies the following recurrence relation 

cosh a; (m — 2) 



Jm ^ X ' = ( m _ 1) sinh™- 1 x (m-1 
which follows from the identity 



(17) 



cosh x 



cosh a; — 



sinh m 2 x ' 



sinh" 1 x sinh m x 

and integration by parts. The antiderivative 3 m (x) naturally breaks into two separate 
classes, namely 



dx 



sinh 2n+1 x 



lV ,+i (2n-l) 
J (2n)!! 



and 



sinh 2n x 



-1 



1 



x . ^ (2A;-2)!!(-l) fc 

log coth — h cosh a; > T , — 

2 (2k- l)!!sinh 2fc x 



^271 -2)!! A (2fe-3)!!(-l) fc 

(2n-l)\\ ^ ^~ nM.^2fc-i 



(18) 



,n+l. 



fc ^ (2fc-2)!!sinh" 
-l) fc coth 2 ^ 1 x 



a 



or 



.r 



(19) 



+ C, 



^ (2k- l)(k- l)!(n- A;)! 
where C is a constant. The double factorial (•)!! : {—1, 0, 1, ...}—> N is defined by 

n • (n — 2) • • • 2 if n even > 2, 
n!! := <( n- (n - 2) • • • 1 if n odd > 1, 
1 if nG {-1,0}. 
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Note that (2n)!! = 2 n n\ for n G N . The finite summation formulae for 3 m (x) all follow 
trivially by induction using ffl7|) and the binomial expansion (cf. (1.2.2) in Olver et al. 
(2010) [23]) 

-l) fc coth 2fc a; 



1 — coth £ 



2 „\n 



?>,! 



E 

fc=0 



fc!(n- k)\ 



The formulae ffl8|) and ffl9|) are essentially equivalent to (2.416.2-3) in Gradshteyn & 
Ryzhik (2007), except (2.416.3) is not defined for the integrand 1/ sinhx. By applying 
the limits of integration from the definition of X^(p) in ( Ti~5"j) to f fTBj) and (jTTJ|) we obtain 
the following finite summation expressions for Z^(p) 

(2fc - 2)!!(-l) fc 



Up) 



" 1 (d-2)\\ 



(_l)(^-l)/2 



log coth — h cosh p 

2 ^ (2£;-l)!!sinh 2fc p 

(d-3)!! 
_(d-2)!! 



if d even, 



d-3 



(d-l)/2 



E 



-lr coth 



2fc-l 



J (2A;-l)(fc-l)!((d-2ife-l)/2)! 



> if d odd. 



or 



/ 1)( d-i)/2 (d-W 
K ' (d-2)\\ 



(d-l)/2 



'V- 2fc-3!!(-l* 

1 + coshp > — i J i 

f-{ (2fc-2)!!sinh 2fe - 1 



(20) 



Moreover, the antiderivative (indefinite integral) can be given in terms of the Gauss 
hypergeometric function 



dp 



-1 



d — 1 d d + 1 



+ C, 



(21) 



sinh^ 1 p (d-1) cosh *" 1 p 2i 1 V 2 ' 2 ' 2 ' cosh 2 

where C G R. The Gauss hypergeometric function 2 -Fi : C 2 x (C \ — N ) x {z G C 
|z| < 1} — )■ C can be defined in terms of the infinite series 

oo 

2 Fi(a,6;c; z) := ^ 



n=0 



{ a )n(b)n z n 

(c) n n! 



(see (15.2.1) in Olver et al. (2010) [2TJ), and elsewhere in z by analytic continuation, 
(see (2.1.5) in Andrews, Askey & Roy 1999), The Pochhammer symbol (rising factorial) 
(•)/ : C -> C is defined by 

n 

(*)n := n^ + z ~ X )' 

i=l 

where / G Nn. Note that 



T(z + l) 
T(z) 



Fundamental solution of Laplace 's equation in hyperbolic geometry 



15 



for all z G C \ — N . The gamma function T : C \ — N — > C (see Chapter 5 in 
Olver et al. (2010) [21]). which is ubiquitous in special function theory, is an important 
combinatoric function which generalizes the factorial function over the natural numbers. 
It is naturally defined over the right-half complex plane through Euler's integral (see 
(5.2.1) in Olver et al. (2010) [21]) 



oo 

2-1. 



T(z) := / t^e-'dt, 
Jo 

Kez > 0. Some properties of the gamma function, which we will find useful are included 
below. An important formula which the gamma function satisfies is the duplication 
formula (i.e., (5.5.5) in Olver et al. (2010) [21]) 

r(2z) = -=-r(z)rlz + -\, (22) 



provided 2z ^ — N , 

The antiderivative ( 12 ip is verified as follows. By using 

-^- 2 F 1 (a, b; c; z) = — 2 F 1 (a + 1, b + 1; c+ 1; z) 

dz c 

(see (15.5.1) in Olver et al. (2010) [21]), and the chain rule, we can show that 

d -1 /d-ldd+1 1 

i Fx 



dp (d - 1) cosh^ 1 p V 2 '2' 2 cosh 2 p 

sinh p p ( ^~ 1 d d+1 1 

cosh d p 2 1 V 2 ' 2 ' 2 ' cosh 2 p 



d sinh p fd+ld + 2d + 3 

+ — r; , 2F1 



(d+l)cosh d+2 p \ 2 ' 2 ' 2 'cosh 2 P/ 

The second hypergeometric function can be simplified using Gauss' relations for 
contiguous hypergeometric functions, namely 

z 2 Fi(a + 1, b + 1; c + 1; z) = -UFAa, b + l\c;z) — iF x (a + 1, b; c; z)] 

a — 

(see p. 58 in Erdelyi et al. (1981) [6]), and 

2 Fi(a, 6 + 1; c; z) = ^-^- 2 F 1 (a, b; c; z) + ^- 2 F 1 (a + 1, b; c; 

(see (15.5.12) in Olver et al. (2010) [21] )• By doing this, the term with the 
hypergeometric function cancels leaving only a term which is proportional to a binomial 
through 

2 F!(a, &;&;*) = (1 - z)~ a 

(see (15.4.6) in Olver et al. (2010) [21]). which reduces to 1/ sinh^" 1 p. By applying the 
limits of integration from the definition of Id{p) m ( TT5T) to ( I2TT) we obtain the following 
Gauss hypergeometric representation 

1 /d- 1 d d+ 1 1 \ 

UP) ~ (d-l)cosl/-V 2Fl l~ 2 ; co^J • (23) 
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Using (1231) . we can write another expression for Id(p)- Applying Eulers's transformation 

2 Fi(a, b; c; z) = (1 - z) c ~ a ~ h 2 F x (c-a,c-b; c; z) 
(see (2.2.7) in Andrews, Askey & Roy (1999) 0), to (122) produces 

MP) - t^^i — ; • , d- 2 „ 2^1 ^ l; -^5—; , 2 



(d — 1) coshp sinh p~ \2 2 cosh p / 
Our derivation for a fundamental solution of Laplace's equation on the i?-radius 
hyperboloid in terms of the associated Legendre function of the second kind is as 
follows. By starting with (1231) and the definition of the associated Legendre function of 
the second kind : (1, 00) — > C, namely 

0Fe^r(^ + //+ l)(z 2 - iy' 2 (v + p + 2 v + p + l_ 3 1 



Q ^ z): 2^(1/ + 1)^+^-1 2Fl { 2 ' 2 ;l/+ 2V 

for |z| > 1 and v + p, + 1 ^ — N (see (8.1.3) in Abramowitz & Stegun (1972) p]), we 
derive 

<m _jn = ^im^v c , (Mhp) . (24) 



2' 2 'cosh 2 py v ^Fe-( rf / 2 " 1 )(d-2)!sinh d/2 - 1 p ^ 



d,0 



We have therefore verified that the harmonics computed in section 13.11 namely u 2 ' + 
give an alternate form of a fundamental solution for the Laplacian on the hyperboloid 
Using the duplication formula for gamma functions ( 1221 . ( 1231) . and ( 1241) . we derive 



e -i7r(d/2-l) 



UP) = 2*-'r(i/2)tfnh*»-' ^K"*")" 



Notice that our chosen fundamental solutions of the Laplacian on the hyperboloid 
have the property that they tend towards zero at infinity (even for the d = 2 case, unlike 
Euclidean fundamental solutions of the Laplacian). Therefore these Green's functions 
are positive (see Grigor'yan (1983) [10]; Grigor'yan (1985) [TT] ) and hence is not 
parabolic. Note that as a result of our proof, we see that the relevant associated Legendre 
functions of the second kind for d G {2, 3, 4, 5, 6, 7} are (cf. (1201 ) 

Qo (coshp) = log coth ^, 

,1/2/ , v . Fk 



sinh 1/2 p ^i/2( cos M = «y 2^ COthp ~ 1 

1^1,,.. , p cosh p 
— (^(cosh p) = log coth — 2 



sinh p 2 sinh p 

-Q3/2 (coshp) = 3zi /— ( — coth 3 p + cothp 



sinh 3 / 2 p^ 3/2V ^ V 2 V 3~"" r r 3 

1 2 , . p cosh p cosh p 

-<y 2 (coshp) = 3 log coth- — 2 : - 4 3 : i 2 , and 



sinh p 2 sinh p sinh p 

-Q5/2 (coshp) = — ( — coth 5 p — - coth 3 p + coth p 



sinh 5 /V 5/2V " ~ V 2 V 15 — ' Z h-t^h 15 
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The constant Cq in a fundamental solution for the Laplace operator on the 
hyperboloid f lT6|) is computed by locally matching up the singularity to a fundamental 
solution for the Laplace operator in Euclidean space, Theorem 13.21 The coefficient cq 
depends on d. It is determined as follows. For d > 3 we take the asymptotic expansion 
for cold(p) as p approaches zero and match this to a fundamental solution of Laplace's 
equation for Euclidean space given in Theorem 13.21 This yields 

r (d/2) 

Co = ~^W' (25) 
For d = 2 we take the asymptotic expansion for 

c X 2 (p) = c logcoth| ~ c log ||x - x'll -1 

as p approaches zero, and match this to {? 2 (x, x') = (2n)~ 1 log ||x — x'|| _1 , therefore 
Co = (2iv)~ 1 . This exactly matches ( 125]) for d = 2. The derivation that Td{p) is an 
fundamental solution of the Laplace operator on the hyperboloid H R and the functions 
for Id(p) are computed above. 

The sectional curvature of a pseudo-sphere of radius R is —1/R 2 . Hence using 
results in Losev (1986) [19], all equivalent expressions in Theorem 13.11 can be used for 
a fundamental solution of the Laplace-Beltrami operator on the i?-radius hyperboloid 

(cf. section namely (where R is now a free parameter) 

The proof of Theorem 13.11 is complete. 

Furthermore, due to a theorem proved in [19], all equivalent expressions for Td(p) in 
Theorem 13.11 represent upper bounds for a fundamental solution of the Laplace-Beltrami 
operator on non-compact Riemannian manifolds with negative sectional curvature not 
exceeding —1/R 2 with R > 0. 

We would also like to mention that a similar computation for a fundamental solution 
of Laplace's equation on the positive-constant sectional curvature compact manifold, the 
i?-radius hypersphere, has recently been computed in Cohl (2011) [5]. 

4. Uniqueness of fundamental solution in terms of decay at infinity 

It is clear that in general a fundamental solution of Laplace's equation in the hyperboloid 
model of hyperbolic geometry 7-L d R is not unique since one can add any harmonic function 
h : — > R to l-L d R and still obtain a solution to 

-A^(x,x / ) = 5 s (x,x / ), 
since h is in the kernel of —A. 
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Proposition 4.1. There exists precisely one C°° -function H : (H R x H R ) \ {(x, x) : 
x G H^,} — > R such that for all x' G tfte function H x / : \ {x'} — > R defined by 



H x i(x) := if(x, x') is a distribution on H R with 



-AH« = 6 g (;x!) 



and 



lim £k(x) = 0, (26) 

a(x,x' J— >oo 

where rf(x, x') zs ifte geodesic distance between two points x, x' G H^. 

Proof. Existence: clear. Uniqueness. Suppose and are two such functions. Let 
x' G H d R . Define the C°°-function ft : H d R \ {x'} ->■ R by ft, = # x / - H x ,. Then ft is 
a distribution on with —Aft = 0. Since H R is locally Euclidean one has by local 
elliptic regularity that ft can be extended to a C^-function ft : H R — > R. It follows 
from (1261) for H and H that 

lim ft(x) = 0. (27) 

d(x,x')— >oo 

The strong elliptic maximum/minimum principle on a Riemannian manifold for a 
bounded domain Q states that if u is harmonic, then the supremum/infimum of u 
in Q coincides with the supremum/infimum of u on the boundary dVL. By using a 
compact exhaustion sequence Qk in a non-compact connected Riemannian manifold 
and passing to a subsequence x^ G dVL^ such that x^. — )• oo, the strong elliptic maxi- 
mum/minimum principle can be extended to non-compact connected Riemannian man- 
ifolds with boundary conditions at infinity (see for instance section 8.3.2 in Grigor'yan 
(2009) [12]). Taking C H R , the strong elliptic maximum/minimum principle for 
non-compact connected Riemannian manifolds implies using fl27j) that ft = 0. Therefore 
ft = and H(x, x') = #(x, x') for all x G \ {x'}. 

By Proposition 14.11 f° r d > 2, the function "H^ is the unique normalized fundamen- 
tal solution of Laplace's equation which satisfies the vanishing decay (1261) . 
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